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Abstract 


A mathematical model was developed which simulates the self-discharge capacity losses in the carbon anode for a SONY 18650 lithium-ion 
battery. The model determines the capacity loss during storage on the basis of a continuous reduction of organic solvent and de-intercalation of 
lithium at the carbon/electrolyte interface. The state of charge, open circuit potential, capacity loss and film resistance on the carbon electrode were 
calculated as a function of storage time using different values of rate constant governing the solvent reduction reaction. 
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1. Introduction 


There are several processes which cause the self-discharge 
capacity loss of lithium-ion batteries during storage, includ- 
ing internal electron leakage, dissolution of active electrode 
materials, corrosion of current collectors and parasitic chem- 
ical/electrochemical reactions on the electrode surfaces. It has 
been known [1-7] that the overall self-discharge of a lithium-ion 
cell is mainly controlled by the parasitic reaction of electrolyte 
reduction or decomposition on both graphite anode and LiCoO2 
cathode. 

The parasitic reactions on the carbon anode have been 
taken into account in the simulation model of charge—discharge 
cycling of lithium-ion cells by Darling and Newman [8]. 
Ramadass et al. [9] and recently Ning and Popov [10] have 
extended the lithium intercalation model to develop a capac- 
ity fade model during cycling process on the basis of the loss 
of active lithium-ions due to the parasitic reactions. Besides the 
fact that there is considerable experimental evidence in the liter- 
ature with regard to the self-discharge, there was no attempt 
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to simulate the self-discharge capacity loss in lithium-ion 
cells. 

All capacity loss of a lithium-ion cell during storage is not 
totally irreversible: namely, a fraction of the capacity lost during 
storage, which is called a ‘reversible capacity loss’ (Qs rey), can 
be recovered upon subsequent charging process, and the remain- 
ing fraction that cannot be recovered is an ‘irreversible capacity 
loss’ (Qsjirr). The underlying mechanism for the ‘reversible 
capacity loss’ is briefly presented in the following section. The 
capacity loss (Qs tot) measured at the end of storage represents 
the sum of Qs rey and Qs irr, and the stored cell should be charged 
and then discharged, in order to experimentally distinguish 
between Qs rey and Qs irr, as explained in the flowchart of Fig. 1. 

The objective of this work is to develop a mathematical model 
which is capable of predicting the total capacity loss occur- 
ring in carbon electrode for lithium-ion cells during storage. No 
attempt was made to distinguish between reversible and irre- 
versible capacity losses, since most applications (e.g. ‘back up 
power systems’ in satellites and space shuttles) may require the 
stored cell to be used ‘as-is’ without re-charging process. 


2. Model development 


The model assumes that the self-discharge reaction in the 
carbon anode involves only the solvent reduction reaction on 
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Nomenclature 

a specific surface area of porous electrode (m~!) 

Chi lithium concentration in the carbon electrode 
(mol m~?) 


Climax Maximum lithium concentration in the carbon 
electrode (mol m7?) 

Cp product (precipitate) concentration on the elec- 

trode surface (mol m~?) 

maximum product (precipitate) concentration on 

the electrode surface (mol m7?) 

Cs solvent concentration at the electrode/electrolyte 
interface (mol m7?) 


CS max Maximum solvent concentration at the elec- 
trode/electrolyte interface (mol m7?) 

F faraday constant (C mol`!) 

Js current density for solvent reduction reaction 
(Am) 

Js0 exchange current density for solvent reduction 
reaction (A m~?) 

ko rate constant for solvent reduction (m s7!) 

Mp molecular weight of newly formed surface film 
(kg mol~ 1) 

Qs irr irreversible capacity loss (C m~?) 

Qsrey reversible capacity loss (C m7?) 

Qstot total capacity loss (C m~?) 

R gas constant (J mol`! K7!) 

Rlm total film resistance (Q m?) 

Rini resistance of initial SEI film (Q m?) 

Rp resistance of newly formed surface film (Q m?) 

t time (s) 

T absolute temperature (K) 

T; total parasitic reaction time (s) 

U local equilibrium potential (V) 

Greek 

a transfer coefficient for solvent reduction reaction 

Ôfilm surface film thickness (m) 

Ns overpotential for solvent reduction reaction (V) 

Kp conductivity of newly formed surface film 
(Sm!) 

pP density of newly formed surface film (kg m~?) 


Subscript or superscript 


a anodic reaction 

c cathodic reaction 

n negative electrode (anode) 
S solvent reduction reaction 


the electrode surface and the resulting loss of active (cyclable) 
lithium-ions. There are no other means of the self-discharge 
such as internal electron leakage, active material dissolution or 
corrosion of the current collectors. 

From the experimental studies on the self-discharge of 
graphite anode, Yazami and Reynier [3,11] have suggested the 
‘metastable adsorbed complex model’ for the self-discharge 


reaction as follows: on the carbon electrode surface, the solid 
phase is in contact with the liquid electrolyte present inside 
the porous solid-electrolyte interphase (SED, as illustrated 
in Fig. 2. During storage of a lithium-ion cell, the elec- 
trolyte percolates through the cracks of the SEI film, and 
lithium-ions tend to de-intercalate slowly from the carbon 
electrode: 


Li? C6 > Liy_y~Co + x(e~ - Ligt) (1) 


where Li®C, is the fully charged carbon anode, Lii—x°°C6 
the partially discharged carbon anode and (e7 -Lig*) repre- 
sents a pair of electron and lithium-ion localized in the 
SEI film resulting from the de-intercalation reaction. Under 
open circuit conditions the lithium de-intercalation is ther- 
modynamically less feasible and is driven solely by the 
reactivity of lithium with the oxidizer present on the carbon 
surface. 

The de-intercalated lithium-ions are consumed by a continu- 
ous formation of metastable electron-ion-electrolyte complex 
during storage. Such a metastable complex dissociates back 
to its components during the charging process of the cell 
after storage, and adsorbed lithium-ions on the graphite edge 
planes are re-intercalated into the carbon electrode. So, this 
reaction is responsible for a ‘reversible capacity loss’ (Qs rev) 
during storage which is recovered upon the charging process. 
On the other hand, a part of metastable complexes would 
undergo the chemical reaction which results in the forma- 
tion of insoluble salts and hence the ‘irreversible capacity 
loss’ (Qs irr). 

As an example, it was assumed in this study that the ethy- 
lene carbonate (EC) is reduced on the carbon anode surface to 
metastable complexes via two-electron transfer during storage 
[12]: 


2EC + 2(e7 - Ligt) > 2(CH2CH2OCODLi*)mmetastable (2) 


Then a part of metastable complexes are transformed into insol- 
uble salts through chemical disproportionation reaction on the 
carbon surface during storage: 


2 y(CH2CH20COzLi” ) metastable 


chemical disproportionation 


== y(CH2)2 + y(CH20C02Li) 4 6) 


The remaining part (1— y) of metastable complexes disso- 
ciate back to the original components upon the following 
charging step of the cell after storage. In this work, the 
model calculates the total capacity loss (Qs,tot) measured 
during storage which is determined by reaction (2). The sol- 
vent reduction reaction is feasible at all temperatures when 
the anode potential is more cathodic than the reversible 
potential for the EC reduction reaction (0.4 V versus Li/Li*) 
[12]. It should be noted that at open circuit potential, no 
current flows through the external circuit and the oxida- 
tion/reduction reactions occur simultaneously within the SEI 
film. 

The following assumptions were made for formulation of the 
model: (i) the solvent reduction reaction is continuous and does 
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Lithium ion cell with an 
initial capacity of ‘Qin’ 


v 
Step 1: 


Lithium-ion cell is stored for 
duration ‘t. 

(Self-discharge process 
occurs during this step) 


R.P. Ramasamy et al. / Journal of Power Sources 166 (2007) 266-272 


Step 2: 


Fully discharge the stored 
cell after time ‘t. 


Step 3: 
Fully re-charge the cell. 


Observation: 


The capacity measured in step 2 will be 
the ‘residual capacity’ (Q_,.) of the cell 
after storage. 


Qies = Qiri - Qs tot 


Observation: 


Mechanism of reversible gain (Q; rev ) 
occurs. A fraction of Q, is recovered 

as Q, ,., during the lithiation (charging) of 
step 3. 


Step 4: 
Fully discharge the cell. 


Observation: 


The actual capacity ( Q,.,) of the cell is 
measured in step 4. Q,., will be less 
than Qin because of the permanent or 
‘irreversible self discharge losses ’. The 
true irreversible capacity loss ( Qù ir ) 
occurred during storage over the period 
of time ‘f is given by 


Qin = Qini = Qact 


Observation: 
Based on the above observations, the 


reversible capacity regained (Q, rev) 
during Step 3 can be calculated as: 


Qs tev = Qstot E Qs irr 


Fig. 1. Flowchart explaining the stepwise procedure for the determination of total, reversible and irreversible capacity losses during storage. 


not depend on the stage transitions in the graphite electrode; 
(ii) the solvent reduction occurs evenly over the entire electrode 
surface; (iii) the solvent is present in excess at the anode/SEI 
film interface and thus does not limit the reduction reaction rate; 
(iv) there is no kinetic limitation for solvent diffusion through 
the cracks of the SEI film; (v) the overall self-discharge rate is 
so slow that it is not limited by lithium diffusion in the carbon 
electrode. 
The rate of the solvent reduction reaction J; on the carbon 
electrode surface can be expressed as 
2 
J = Ja l os 
n) | 


2 
Cs „max Cii, max 


exp ( 


Cp 
Cp, max 


Qa F 
RT 


(4) 


According to the assumption (v), lithium is uniformly distributed 
throughout the carbon particle: namely, the surface concentra- 
tion of lithium is equal to the bulk concentration inside the carbon 
electrode. In Eq. (4), the exchange current density for solvent 
reduction Jso depends on the initial lithium concentration in the 
carbon electrode: 
Jso = nFkgC{ (5) 
Here, the rate constant for solvent reduction ko is dependent 
on the electrolyte composition and electrode structures in the 
specific battery as well as the storage temperature. The overpo- 
tential for solvent reduction ns under open circuit conditions is 
expressed as: 

= USC _ yocr 


Ns (6) 
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Surface Film (SEI) 


T 


Solvent reacts with the de-intercalated 
lithium ion and electron and forms product 


& Solvent molecule 
E Product 


G Solution phase lithium ion 


De-intercalation 
of lithium at the 
surface 


Fig. 2. Schematic diagram of the self-discharge mechanism on the carbon electrode in a lithium-ion cell. 


The discharge curve of the carbon electrode taken from a 
SONY 18650 lithium-ion cell was measured at a C/50 rate 
in a three-electrode electrochemical cell with the lithium 
reference and counter electrodes, and then the functional 
expression of uo was obtained by a non-linear fit- 
ting of the experimental discharge curve, as given in the 
Appendix. 

Since the potential window of interest is far from the 
reversible potential for EC reduction, it is reasonable to model 
the solvent reduction reaction with Tafel kinetics [8,9] as fol- 


lows: 
J J Cs i ( æF ) (7) 
= Js0 exp n 
i i Cs, max Cin RT j 
Based on the assumptions (iii) and (iv), Eq. (7) reduces to 
CLi 4 ac F 
Jg=Jy —— } exp | —-——», 8 
s s0 5) p RT Ns (8) 


Since the lithium concentration depletes as fast as the solvent 
reduces, the following boundary condition was used for numer- 
ical calculation: 


(9) 


The loss of active lithium-ions Qs,tot per unit surface area of 
carbon particle during storage was estimated using the following 
equation: 


Ts 
Qs tot = | Js dt 
0 


The total resistance of the SEI film, Rfim, on the carbon electrode 
is the sum of the time-invariant resistance Rini of the initial SEI 
film formed during the formation period and the time-dependent 
resistance Rp of the SEI film which is newly formed during 


(10) 


Table 1 
Key parameters used in the simulation model [9] 


Parameter Unit Value 


m7! 


Specific surface area of porous 3.0 x 10° 
electrode, a 

Maximum lithium concentration 
in the carbon electrode, Climax 

Molecular weight of surface film, 
Mp 

Reversible potential for solvent 
reduction, UOCP 

Transfer coefficient for solvent - 
reduction, œe 

Conductivity of surface film, xp 


Density of surface film, pp 


mol m~? 3.056 x 10* 


kg mol7! 7.3 x 10+ 
V vs. Li/Li* 0.4 


0.5 


storage: 


Rfm(t) = Rini + Rp(t) (11) 


where 


1 of 
Rp(t) = a 


(12) 


The thickness of the SEI film increases continuously over time 
due to the electrolyte reduction reaction: 


OS fim = J;Mp 
ot nopF 


(13) 
The key parameters used in the model are listed in Table 1. 
3. Results and discussion 

The model predictions were performed based on the experi- 


mental data obtained for the carbon anode taken from a SONY 
18650 lithium-ion cell [13,14]. The electrochemical exper- 
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Fig. 3. Plots of OCP of the carbon electrode against storage time measured 
experimentally (symbol) and simulated theoretically (line). The simulations 
were performed with kọ =1.5 x 107!8 and 1 x 107!7 ms! at 25 and 35°C, 
respectively. 


iments were conducted in a three-electrode electrochemical 
cell (T-cell) with the lithium reference and counter electrodes. 
The electrolyte used was 1 M LiPFe-EC/PC (50:50 vol%). Ini- 
tially, the carbon electrodes were fully lithiated at 30mV 
versus Li/Lit, followed by 1-h rest to ensure that the potential 
reaches its equilibrium value. Next, the OCPs were monitored 
as a function of storage time at different temperatures of 25 
and 35°C. 

Fig. 3 shows the plots of the open circuit potential (OCP) 
of the carbon electrode against storage time measured exper- 
imentally and simulated theoretically. In the model, the rate 
of solvent reduction reaction Js on the carbon electrode was 


(a) 100 


SOC of Carbon / % 


20 


k,=1x10" ms" an 


Storage Time / month 
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assumed to be given by the Tafel equation (Eq. (8)), and the 
value of rate constant kọ was used as an adjustable parameter. 
As shown in Fig. 3, it was found that the OCP values calculated 
with kọ = 1.5 x 107!8 and 1 x 10717 m s7! give the best fit to the 
experimental data for 25 and 35 °C, respectively. In fact, such ko 
values yield the initial exchange current densities Jso|;-9 which 
are in agreement with those given in literature [9]. The model 
does not fit the experimental data exactly at time t=0. The ini- 
tial condition used for the model is 100% state of charge (SOC) 
for the carbon electrode. In practice, however, the carbon elec- 
trode cannot be fully lithiated to its theoretical maximum value, 
and hence the experimental data at t=0 show higher OCPs than 
those predicted by the model. 

The lithium concentration in the carbon electrode was deter- 
mined as a function of storage time by running the model with 
various values of ko, and the simulation results were presented 
in Fig. 4(a) for a period of 1 year. The carbon electrode, before 
undergoing the self-discharge, was completely lithiated, i.e. 
SOC = 100%. As shown in Fig. 4(a), the self-discharge rate of 
the carbon electrode is very high in the initial stage of storage 
and then levels off for longer storage time. The SOC decreases 
with storage time faster for higher ko value, due to higher reactiv- 
ity of lithium with the solvents. Under open circuit conditions, 
the potential of the carbon electrode is determined solely by 
the solid phase concentration of lithium. Fig. 4(b) shows the 
OCP values of the carbon electrode simulated as a function of 
storage time for different kọ values. The OCP value increases 
more rapidly with storage time for higher kọ value due to higher 
de-intercalation rate of lithium from the carbon electrode, as 
described in Fig. 4(a). 

Fig. 5 illustrates the OCP of the carbon electrode simulated 
as a function of SOC for 5 years of storage. The simulation 
predicts a decrease of SOC from 100 to 9% and an increase 
of OCP from 0.047 to 0.249 V versus Li/Li* after 5 years of 
storage. 


(b) 


OCP of Carbon / V vs. Li/Li” 
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Storage Time / month 


Fig. 4. Plots of (a) SOC and (b) OCP of the carbon electrode against storage time simulated with different kg values. 
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Fig. 5. OCP of the carbon electrode simulated as a function of SOC for 5 years 


of storage. The simulation was performed with kọ =1 x 107!7 ms7!. 


Fig. 6 presents the solvent reduction rate J; calculated as a 
function of storage time using different kg values. As expected, 
the higher the ko value is, the higher is the initial J; value; how- 
ever, J; decreases with storage time at a faster rate for higher ko 
value. For storage periods longer than 15 days, the J, value for 
ky =1 x 107!8 ms7! is lower than that for kg = 1 x 1071? ms7!. 
This is due to the fact that for higher kọ the lithium concentration 
in the carbon/electrolyte interface depletes at a faster rate and 
reaches very low values within a short period of storage time, as 
shown in Fig. 4(a). 

Fig. 7 shows the capacity loss Qs,tot and the film resis- 
tance Rp simulated for a period of 1 year as a function 
of storage time using different kọ values. The model pre- 
dicts a continuous increase in the capacity loss and film 
resistance with storage time. From Eqs. (10) to (13), it is 


N -1 
g 10 
(5) 
< 
E 10? 
% 
oO 
w 10° 
X 
© 
2 1 : 
Ə 10* $ 0.05 0.10 0.15 0.20 
oD — k =1x10" ms" 
w 51 _ -18 4 
€ 10 K k =1x10 ms 
v - - 
2 k,=1x10" ms" 
D 498 

107 

0 2 4 6 8 10 12 


Storage Time / month 


Fig. 6. Plots of solvent reduction rate J; against storage time simulated using 
various kg values. 
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Fig. 7. Plots of capacity loss Qs,tot and film resistance Rp against storage time 
simulated for different kọ values. 
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Fig. 8. Plots of (a) capacity loss due to the self-discharge of the carbon elec- 
trode and (b) overpotential for solvent reduction as a function of storage time, 
simulated for different initial SOC values. 
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noted that Rp is linearly proportional to Q¢ tor in the carbon 
electrode. 

Fig. 8(a) presents the capacity loss Qs to due to the self- 
discharge of the carbon electrode simulated for different initial 
SOCs. For the same value of kọ =1 x 107!” ms7!, the capacity 
loss increases with increasing the initial SOC value. The results 
can be explained by taking into account that the OCP of the car- 
bon electrode UOC? with higher SOC is more cathodic to the 
solvent reduction potential Ue than that of the electrode with 
lower SOC. In other words, the overpotential for the solvent 
parasitic reaction increases with an increase of the SOC of the 
carbon electrode, as presented in Fig. 8(b). 


4. Conclusion 


The mathematical model developed in this study predicts the 
capacity loss of the carbon electrode during storage due to a 
continuous reduction of organic solvent and de-intercalation 
of lithium at the carbon/electrolyte interface. The SOC, OCP, 
capacity loss and film resistance on the carbon electrode were 
simulated as a function of storage time using different values 
of rate constant for solvent reduction kp. The effect of initial 
SOC on the capacity loss was also studied. The simulated values 
of OCP were found to be in good agreement with the exper- 
imental values obtained for the carbon electrode in a SONY 
18650 lithium-ion cell. In the following work, the self-discharge 
model presented in this study will be combined with the capacity 
fade cycling model [10], in order to analyze the calendar life of 
lithium-ion batteries during intermittent cycling process. 


Appendix A 


The UOC expression for the carbon electrode taken from a 
SONY 18650 cell is given by 


OCP CLi Cu \" 
U; = 0.7222 + 0.1387 | ——— ] + 0.029| ——— 
Li,max CLi,max 


Cr; =l Ci: —1.5 
=) + o.019( =) 
CLi,max CLi,max 


C 3 
+ 0.2802 exp fos —15 (+) \ 
CLi,max 


Chi 


-00172( 


— 0.7984 exp {0.4405 ( ) — 04108} (A.1) 


Li,max 
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